Nontrivial critical behavior of the free energy in the two-dimensional Ising spin glass 

with bimodal interactions 
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A detailed analysis of Monte Carlo data on the two-dimensional Ising spin glass with bimodal 
interactions shows that the free energy of the model has a nontrivial scaling. In particular, we 
show by studying the correlation length that much larger system sizes and lower temperatures are 
required to see the true critical behavior of the model in the thermodynamic limit. Our results agree 
with data by Lukic et al. in that the degenerate ground state is separated from the first excited 
state by an energy gap of 2 J. 
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I. INTRODUCTION 

The two-dimensional Edwards- Anderson Ising spin 
glass with bimodal interactions! has been widely used to 
study properties of spin-glass systems. Despite the fact 
that the model only orders at zero temperature, its pop- 
ularity can be ascribed mainly to its ease of implemen- 
tation and simplicity. There has been a long-standing 
controversy regarding the behavior of the free energy of 
the model, in particular the size of the excitation gap and 
the critical behavior in the thermodynamic limit. 

The exponential scaling of the free energy, and thus 
correspondingly of all other thermodynamic quantities, 
was first proposed by Wang and Swendsen. 2 They sur- 
mised that 
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where j3 = 1/T represents the inverse temperature, J is 
the magnitude of the bonds, and A is a numerical pref ac- 
tor. Simple analytical arguments for Ising systems with 
bimodal interactions on square lattices with coordination 
4 suggest that the energy gap should be 4 J, i.e., A = 4. 
In addition, according to hyperscaling, the singular part 
of the free energy scales as with d = 2 and so, if 
Cy scales exponentially, we expect that the correlation 
length £ scales as 



(2) 



with A = 2n, as predicted first by Saul and Kardar^ 

Wang and Swendsen were the first to calculate numer- 
ically the specific heat of the model and find A = 2, thus 
showing a nontrivial scaling of the free energy; however 
their results for small system sizes and few disorder re- 
alizations implied strong corrections to scaling. These 
results were later supported by numerical work of Lukic 
et aiA who computed the specific heat of the model for 
intermediate system sizes (L < 50) and also showed that 
A = 2. 



The aforementioned results stand in contrast to work 
by Saul and Kardar^who argue that A = 4. In addition, 
in their original work n = 2, a behavior later confirmed 
independently by Houdaycr— who studied the finite-size 
scaling of the Binder ratio^ as well as Katzgraber and 
Leei who studied the finite-size scaling of the finite-size 
correlation length directly via Monte Carlo simulations. 
Note that exponential scaling of the correlation length 
means that the critical exponent v — oo for the two- 
dimensional spin glass with bimodal interactions. 

In this work, by using an alternate analysis of the data 
presented in Ref . for the finite-size correlation length^ 
we show that n changes continuously with system size. In 
particular, our analysis shows that n — 1 cannot be ruled 
out (in contrast to previous conclusions - see Ref. @) and 
might be the proper estimate for n in the thermodynamic 
limit. In addition, we compute the specific heat for mod- 
erate system sizes directly via Monte Carlo simulations 
and show that A — 2, in agreement with data by Lukic 
et al£ 

In Sec. [H] we introduce the model and observables and 
in Sec. IIHI we present our results, followed by conclusions 
in Sec. US 



II. MODEL AND OBSERVABLES 

The Hamiltonian of the two-dimensional Ising spin 
glass is given by 
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Si = ±1 represent Ising spins and the sum is over near- 
est neighbors on a square lattice with periodic bound- 
ary conditions. The interactions 6 {±1} are bi- 
modally distributed. For the Monte Carlo simulations 
we use a combination of single-spin flips, parallel tem- 
pering updatesfi2iii and rejection-free cluster moves^ to 
speed up equilibration. Equilibration of the method is 
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TABLE I: Parameters of the simulations. AT samp represents 



the number of disorder realizations computed; JV S , 



is the 



total number of Monte Carlo sweeps of the 2Nt replicas for a 
single sample. Nt is the number of temperatures in the par- 
allel tempering method and T m in represents the lowest tem- 
perature simulated. 
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N', sweep 




A?T 


32 


5 000 


2.0 x 10 6 


0.050 


20 


48 


1 000 


2.0 x 10 6 


0.050 


20 


64 


500 


4.2 x 10 6 


0.200 


39 


96 


609 


6.5 x 10 6 


0.200 


63 


128 


420 


2.0 x 10 6 


0.396 


50 



tested by performing a logarithmic data binning of all 
observables, and we require that the last three bins agree 
within error bars and are independent of the number of 
Monte Carlo sweeps N swee p. The parameters of the sim- 
ulation are listed in Table ||] 

The finite-size correlation length 8 i 9 i 12 i 13 i 14 i 15 £l is 
given by 
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where k m j n = (2ir/L, 0) is the smallest nonzero wave vec- 
tor, and xsc(k) is the wave- vector-dependent spin-glass 
susceptibility: 



XsG (k) = l^[(5^ 3 ) 2 ] av e^( R -^ 



(5) 



In the previous equation [• • -] av represents a disorder av- 
erage and (• ■ ■) a thermal average. 

In order to compare to the data of Lukic et al£, we 
also compute the specific heat of the system: 
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III. RESULTS 



In Fig. ^ w e show data for the natural logarithm of 
the finite-size correlation length as a function of 1/T for 
different system sizes. We expect data for 1ii(£l) vs 1/T 
to approach a slope of n. The data, at first sight, show 
good agreement with n — 2 (dashed line). Note that the 
data peel off the asymptotic behavior thus masking any 
potential effective variation of n with L and T. Data at 
high temperatures show a slope n > 2, thus suggesting 
that the slope might change for decreasing temperatures 
and increasing system sizes. 

Figure [2] shows data for the derivative of the natural 
logarithm of the finite-size correlation length with respect 
to the inverse temperature as a function of temperature 
for several system sizes. In the bulk regime (T > 0.7) 
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FIG. 1: (Color online) Natural logarithm of the finite-size 
correlation length as a function of 1/T for several system sizes 
(data taken from Ref. Because ~ exp(n/3J) we expect 
the slope of the curves in the plot to determine n. For large 
T the data are independent of L and "peel off" from this 
common curve at a value of T which decreases with increasing 
size. The slope of the common curve asymptotically seems to 
approach the value n — 2 (dashed line). The dotted line has 
slope n — 1 which corresponds to the nontrivial scaling. The 
data seem incompatible with this prediction. 



where the data show no system size dependence, the com- 
puted values for the slope suggest 2 < n < 3. Interest- 
ingly, the data vary linearly with T and then successively 
peel off the linear behavior at increasingly lower temper- 
atures for increasing system sizes. A linear extrapolation 
of the data in the bulk regime seems plausible and shows 
that 7i = 1 at low temperatures in the thermodynamic 
limit, thus concluding in a nontrivial scaling of the free 
energy in the two-dimensional Ising spin glass with bi- 
modal interactions. 

Our results clearly show that the two-dimensional Ising 
spin glass with bimodal interactions has strong correc- 
tions to scaling. In Ref. a finite-size scaling analysis 
of the finite-size correlation length suggested compatibil- 
ity with n = 2, whereas a finite-size scaling of the data 
for n = 1 did not seem plausible. Merely the data for 
the two largest sizes showed a signature of a scaling be- 
havior when n = l. This in turn lead to the conclusion 
that n = 2. While the same observable is being studied 
here, the results presented contradict these findings and 
show that to understand the low-temperature properties 
of this model, the limit of L — > oo has to be taken before 
extrapolating to T = 0. 

Finally, in Fig. [21 we show data for the specific heat Cy 
as a function of temperature plotted as — T\n(T 2 Cy) vs 
T. The data plotted in this way should tend to A for 
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FIG. 2: (Color online) Derivative of the natural logarithm 
of the finite-size correlation length with respect to the inverse 
temperature as a function of temperature for several system 
sizes. In agreement with Fig.^the slope for the different data 
sets ranges between 2 and almost 3 for the highest tempera- 
tures, suggesting that the slope n is temperature dependent. 
For large temperatures the data are independent of L and 
well fitted by a straight line (solid line in plot), thus suggest- 
ing that corrections to n are quadratic. In addition, the data 
peel off at low temperatures from the straight-line behavior. 
An extrapolation of the data to T — agrees well with n = 1 
in the thermodynamic limit. Note that the data show that 
the results in the thermodynamic limit are incompatible with 
n = 2 for L — » oo. The dotted lines correspond to error 
estimates to n. 

T — > 0. Our results show that A « 2 for T — > already 
for intermediate system sizes, in agreement with results 
by Lukic et alA 

IV. CONCLUSIONS 

To conclude, we have shown data for the correlation 
length and the specific heat of the two-dimensional Ising 
spin glass with bimodal interactions. The data for the 
correlation length scale exponentially ~ exp(n/3J) and 
seem compatible with n = 1 in the thermodynamic limit, 
in contrast to previous calculations £ A detailed study of 
n as a function of temperature shows that n strongly 



depends on system size and temperature. This result is 
supported by data for the specific heat which also scale 
exponentially with A ks 2 = 2n. This means that the 
excitation gap for the bimodal spin glass is « 2 J. Our 
results in the thermodynamic limit agree with work by 
Lukic et al£ as well as Wang and SwendseniS 

Our data show that corrections to scaling in the two- 
dimensional Ising spin glass with bimodal interactions are 
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FIG. 3: (Color online) Specific heat Cv plotted as 
— Tln(r 2 Cv) vs T for intermediate system sizes. For 
large enough system sizes L and low enough T, data for 
— Tln(r 2 Cv) should tend to A. In this case the data agree 
with the results of Lukic et all— in that A ~ 2 and not A ~ 4. 

extremely large. In particular, it is important to take the 
limit L — > co before extrapolating to T = 0. Thus system 
sizes of at least ~ 100 2 and larger are required to probe 
the true thermodynamic behavior. Therefore an analy- 
sis with yet larger system sizes at lower temperatures is 
imperative. 
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